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ABSTRACT 


Algorithms for blind equalization and data recovery of orthogonal frequency-division 
multiplexed (OFDM) signals transmitted through fading channels are implemented and 
simulated in this thesis. The channel is estimated without knowledge of the transmitted 
sequence (i.e., blindly) using a least mean squares (LMS) adaptive filter and filter bank 
precoders. This method was used to estimate channel characteristics using both binary 
and quadrature phase-shift keying signals. Additionally, the method was analyzed for 
robustness with a poor initial estimate of channel characteristics, with the addition of 


white Gaussian noise to the signal, and with non-stationary channel conditions. 


In addition, it was shown that the proposed method is particularly suited in 
situations with deep fading channels, where some of the subcarriers have a very low 


signal-to-noise ratio (SNR). 


Simulations for both aspects of this thesis were conducted using MATLAB, and 


the results are presented. 
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EXECUTIVE SUMMARY 


Orthogonal frequency-division multiplexing (OFDM) modulation has become the 
preferred modulation method for numerous modern standards, from IEEE 802.1 1a/g/n 
(Wi-Fi) to IEEE 802.16 (WiMAX) to 4G Long Term Evolution [1]. The reason that 
OFDM modulation has become so prevalent is because of a number of advantages it has 


over other modulation techniques. 


One of the main advantages of OFDM modulation is the ease of channel 
equalization, especially in multipath environments [1]. This is achieved by dividing the 
channel bandwidth into narrow sub-channels with orthogonal subcarriers, and 
transmitting a M -ary quadrature amplitude modulation (MQAM) symbol on each of the 
subcarriers. Even in the presence of severe multipath, each subchannel is affected by 
flat-fading, which can be easily equalized. A guard time interval between symbols 
prevents inter-symbol interference. Additionally, OFDM modulation supports high peak 


data rates and non-line-of-sight applications [1]. 


As OFDM symbols are transmitted through a channel, they undergo distortion due 
to noise and multipath effects. To counteract the effects of channel distortion, an 
estimate of channel characteristics must be calculated. In most standards, this is done 
using preambles and pilot signals, resulting in a loss of data rate. Additionally, most 
standards use a cyclic prefix (CP) during their guard intervals, which results in a loss of 


power efficiency since this information is disregarded at the receiver. 


A method of blind channel equalization that can be applied to OFDM has been 
proposed by P. P. Viadyanathan et al [2], [3]. This method uses a zero prefix (ZP) vice a 
CP, the reason being that the transmitted ZP constrains the received signal to a subspace 
from which the channel information is extracted. This increases both data throughput and 


power efficiency. 


The method of blind equalization proposed in [2] is investigated in this thesis. 


The method is implemented in MATLAB with both off-line and real-time channel 


XV 


estimation investigated. Although the method in [2] uses a ZP, we show that it can be 
extended to the most traditional CP based OFDM by some simple preprocessing of the 


received data. 


The method proposed in [2] was implemented as a recursive algorithm using a 
least mean squares (LMS) adaptive estimation method. However, it has been found that 


it is very sensitive to the SNR of the received signal. 


In this research, the approach in [2] has been re-derived and a third, alternative 
adaptive algorithm implemented. This method also uses a LMS adaptive estimation 
method; however, it is shown to be more robust. Both adaptive methods can be fine- 
tuned by properly setting their step size. This thesis investigated two LMS step sizes, 


determined empirically. 


The performance of these algorithms under different conditions is also compared. 
The first simulation was a comparison of their performance in a channel with additive 
white Gaussian noise (AWGN). The third method of channel estimation with the larger 
of the two step sizes performed the best in high signal-to-noise ratio (SNR) conditions 
(around 30dB). The second simulation used a non-stationary channel without AWGN 
modeled as a simple random walk. In this case, the original method of channel 
estimation and the alternative adaptive method had similar results, although the 
alternative adaptive method with the larger step size performed slightly better at channel 


rates of change of 0.08 percent of the total channel power. 


Since both adaptive methods require an estimate of the channel as an initial 
condition, the final test run on the two algorithms was a comparison of how long it would 
take the channel estimate to converge with an inaccurate initial channel estimate. In this 
case, the more robust alternative adaptive method with the larger step size once again 


outperformed the initial adaptive method. 


Finally, this thesis investigates a method of estimating data lost in nulls caused by 
deep fading in the wireless channel, due to multipath. The algorithm was implemented in 


MATLAB, and compared to a normal OFDM receiver that does not conduct null 


XVi 


estimation. The proposed null-estimating algorithm has similar error rates to the original 
receiver for SNRs less than approximately 21 dB but then significantly outperforms it 


and recovers the data lost in nulls for larger SNRs. 


[1] J. G. Andrews et al., Fundamentals of WiMAX. San Francisco, CA: Prentice Hall, 
2007. 


[2] P. P. Vaidyanathan and B. Vrcelj, “A frequency domain approach for blind 
identification with filter bank precoders,” in Proc. 2004 Int. Symp. Circuits and 
Systems, pp. HI- 349-52, 2004. 


[3] P. P. Vaidyanathan and B. Su, “Remarks on certain new methods for blind 
identification of FIR channels,” in Conf: Rec. 38" Asilomar Conf. Signals, 
Systems and Computers, Asilomar, CA, pp. I- 832—836, 2004. 
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I. INTRODUCTION 


A. BACKGROUND 


Orthogonal frequency-division multiplexing (OFDM) modulation has become the 
modulation method of choice for wireless broadband communications. It is currently in 
use in numerous standards, including IEEE standards 802.11 and 802.16, and military 
wireless communications such as the Wideband Networking Waveform of the Joint 
Tactical Radio System. The reason for its proliferation comes from the advantages that it 
has over other methods of modulation. The main advantage is that the technique allows 


high data rate transmissions to be sent over channels that experience multipath fading. 


One of the challenges of this modulation technique is that it requires knowledge 
of the multipath channel to perform equalization. How to acquire knowledge of the 
channel using blind equalization is examined in this thesis. This is significant since it 
increases the throughput of the channel by reducing (or eliminating) the need for pilot 


signals normally embedded in the OFDM symbol. 


B. OBJECTIVES AND METHODOLOGY 


The goal of this thesis was to further investigate the method of blind channel 
equalization proposed by P. P. Vaidyanathan, et al. [1]. Once the method had been 
implemented in code, a binary or quadrature phase-shift keyed (BPSK or QPSK 
respectively) signal was sent through either a channel with additive white Gaussian noise 


(AWGN) or a non-stationary channel modeled by a simple random walk. 


Once this was completed, the method was implemented using a least mean 
squares (LMS) estimation technique. This method was also tested using BPSK and 
QPSK signals through either a channel with AWGN or a non-stationary channel. It was 
also tested using an inaccurate initial channel estimate to determine how long the channel 


estimate would take to converge to the correct channel. 


Additionally, a secondary LMS method was implemented to overcome ill- 
conditioning issues with the initial LMS method. This was tested in the same manner as 


the first LMS method. 


Finally, this paper shows an additional advantage of the secondary LMS method; 
namely, that this technique allows for the recovery of information sent through deep 


fading channels. 


Cy BENEFITS OF STUDY 


Demand for data, specifically multimedia content, continues to increase in a 
nearly exponential manner [2]. This has led to an increased demand on mobile 
broadband technologies, including 4G phone technology and WiMAX [3]. Both these 
emerging technologies use OFDM modulation. By reducing or eliminating pilot signals 
in the OFDM symbol, the amount of data throughput is increased. Additionally, if a zero 
prefix (ZP) is used instead of a cyclic prefix (CP), it saves power at the transmitting 


station, whether it is a mobile or base station [3]. 


D. ORGANIZATION 


This thesis is organized into six chapters. Chapter II is a background discussion 
of OFDM modulation, and examples of its use in current standards are provided. Chapter 
If contains an in-depth discussion of the blind equalization method proposed in [1]. An 
adaptive estimation method based on the discrete Fourier transform (DFT) is proposed in 
Chapter IV. Simulation results from the adaptive method are compared with simulation 
results from the original in Chapter TV. A method of recovering the information lost in 
deep fading channels is discussed in Chapter V, and the results of the implementation of 
this method are presented. Finally, a summary of the work in this thesis, significant 


results, and recommendations for future work are given in Chapter VI. 


HW. OVERVIEW OF OFDM 


Orthogonal frequency-division multiplexing is a method of multicarrier 
modulation used in a variety of standards for both wired and wireless communications. 
A brief overview of OFDM modulation and its use in IEEE standards 802.11a, 802.11g, 
and 802.11n is provided in this chapter. 


The principle behind OFDM is the idea that a wideband frequency channel can be 
divided into subchannels with narrower bands on which data can be transmitted. 
Orthogonality of the subcarriers, in a sense specified later, prevents interference between 
the subbands. This technique is particularly useful when transmitting through frequency- 
selective channels arising from environments that have severe multipath, and it allows 


simple channel equalization at the receiver [4]. 


In OFDM, the data is transmitted in blocks via a sequence of OFDM symbols, 
separated by a guard time interval to prevent inter-symbol interference (ISI). This is 
guaranteed by choosing the guard time between OFDM symbols to be larger than the 


maximum channel time spread [4]. 


A key issue addressed in this thesis is what is best to transmit during the guard 
time. In the most common implementations of OFDM, such as in IEEE Standards 
802.11, 802.16, and others, the CP is used, which is a cyclic repetition transmitted as the 
guard interval. A consequence of choosing to transmit the CP is that it becomes very 
easy to demodulate the signal in the frequency domain, since linear convolution becomes 


circular convolution where the fast Fourier transform (FFT) plays a dominant role. 


In a number of papers [1][5], it is shown that a ZP, obtained by transmitting all 
zeros during the guard interval, is an attractive alternative to a CP because it yields 
sufficient information to perform blind equalization. In other words, in principle at least, 
with the ZP the channel can be estimated from the received data without the need of 
pilots and/or training signals. This is due to the fact that, with the addition of the 
transmitted zeros, the received signal is constrained on a subspace that yields the channel 


information, regardless of the transmitted data. 
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This research takes this issue even further. It shows that another advantage of 
using an OFDM symbol with ZP (rather than the standard CP) is the fact that by using a 
longer FFT at the receiver, subcarriers with low signal-to-noise ratio (SNR) due, for 
example, to deep fading, can still be reliably recovered using simple properties of the 
FFT. Of course, this comes with some added complexity but is still manageable with 


today’s technology. 

Interestingly, with some minor preprocessing, the ZP OFDM approach can be 
extended to the existing CP OFDM with the same advantages of blind equalization and 
robust recovery in deep fading channels. 
as OFDM SYMBOL 


In OFDM, the transmitted signal is defined as 


s(t) =Re{e?"*'x(y} (1) 
where F’, is the carrier frequency, and the baseband complex signal is defined as 
x(t) = >: geen 0 Sts dare (2) 
Ne 


where the coefficients c, depend on the transmitted data and N,, is the number of 


subcarriers. The symbol length T 


symbol 


is made of the data length 7, and the time guard 
length 7, as 


T. 


symbol al 


se ae (3) 
To ensure that the subcarriers are orthogonal to one another, the difference 


between subcarrier frequencies is chosen to be 


eae (4) 


This results in the following integral, which defines the orthogonality of the 


subcarriers: 
1 an j27F, j27F, 1 me j22(k-l)AF 1 ifk =I 
— [ete Paridg=— [| ePeb nrg) * (5) 
T, i, foe 0 ifke#l 
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If the channel is linear time-invariant (LT), at least within the symbol duration, 


and the guard time 7, is greater than the channel time spread, the subcarriers are still 
orthogonal at the receiver and are given by 


y(t)= SGiH(h je T Sted., 
k 


symbol ? 


(6) 
where H(F) is the frequency response of the channel [4]. 


The discrete time implementation is very straightforward. Let F be the sampling 


frequency, N be the number of data samples in each symbol, and AF =1/NT, = F./N be 
the subcarrier spacing. Then (2) becomes 


Nr vs Ne 
1 2 j2ak—(n-L) 1 2 ent) 
x(nT,) = — sy ce % =— Y ce.” n=0,....L+N-1 
N ie NE N pee 
2 


(7) 
2 
where the guard time is 7, = LT, and the scale factor of 1/N is immaterial but is included 


for later convenience. This result leads itself to using the inverse fast Fourier transform 
(IFFT) to complete the precoding necessary for OFDM: 


(8) 
If X[k]=c, 


for positive subcarriers (k>0), 


X[N+k]=c, for negative 
subcarriers (k <0), and X[k]=0 for k =O, then (8) can be written in terms of the IFFT 
as 


x[n+L]= > xtkle = IFFT {X{k]}. 


(9) 


From (2) and the subcarrier frequencies being multiples of AF in (4), it is easy to 
see that the data transmitted during the guard interval is just a periodic repetition. In 
particular, the CP is comprised of the last L points of the FFT {X [k]} : 


The demodulation of the OFDM signal is based on the convolution of the channel 
impulse response and the signal: 


ae, (10) 
Be A[k|X[k]e % 
= IFFT { H[{k]X[k]} 
Based on the results of (10), the following conclusion can be drawn: 
H[k]X[k] = FFT { y[L],..., y[L+ N -1]} | ai 
with H[k]= FFT {h{0],...,h[L-1],0,...,0} k=0,...,N—-1 


The block diagram of a simple OFDM communication system is shown in Figure 
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Figure 1. A block diagram of a simple OFDM communications system. From [4]. 


Since for the k-th subcarrier, the received signal is of the form 
Y[k]= A[k|X[k]+W[k] (12) 
where W[k] is the FFT of win], additive noise, the transmitted signal X[k] can be 


recovered using a Weiner Filter as 


H'[k] 


X[k] =——— 
|H[k]| +o 


(13) 


if the channel frequency response is known. 


There are some limitations to this approach of OFDM modulation. The CP does 
not contain useful data and is just overhead for the modulation technique. Additionally, 
the channel imposes some constraints on the length of the data interval such that it must 
be larger than the maximum channel time spread to minimize overhead and must be 
much smaller than then period associated with the Doppler spread to maintain the LTI 


channel assumption [4]. 


The advantages of OFDM are that it is effective in multipath environments, it is 
efficient in implementation, and the data rate can be optimized based on channel 
conditions and noise characteristics. The disadvantages are that it is sensitive to time and 
frequency synchronization, and data can be lost when there are spectral nulls in the 


channel frequency response [4]. 


B. OFDM IN STANDARDS 
1. IEEE Standard 802.11a 


TEEE Standard 802.1la, which describes the high-speed physical layer in the 
5GHz band, uses OFDM to transmit data over wireless local area networks (LANs) [6]. 
The standard uses 52 sub-carriers, four of which carry the pilot signal and 48 that carry 
the data. The AF is set at 0.3125 MHz, with the guard interval 0.8 ws long and the data 
interval 3.2 ws long. The mapping of the 52 sub-carriers to the time domain outputs can 


be seen in Figure 2. 


Null 


Time Domain Outputs 





Figure 2. Mapping of IFFT inputs to time-domain outputs for 802.1 1a. From [7]. 


2. IEEE Standard 802.11g 
IEEE Standard 802.11g, is a modification of the 802.11 standards for the 2.4-GHz 


band to include OFDM modulation. The inputs to the IFFT retain the same structure as 


those in 802.11a but with a lower carrier frequency 


3. IEEE Standard 802.11n 


TEEE Standard 802.11n uses OFDM modulation in a multiple-input multiple- 
output (MIMO) scheme to attain higher throughput and more robustness to multipath and 
noise for both the 2.4- and 5-GHz bands [8]. It provides up to four spatial streams of 


OFDM modulated data to achieve a maximum data rate of 600 Mbps. 


4. Other Standards 


OFDM modulation is also used in various other applications. WiMAX is based 
on IEEE Standard 802.16, an OFDM based system. Other standards that use OFDM 
modulation include power-line communications, various digital television and radio 


standards, and the Wideband Networking Waveform of the Joint Tactical Radio System. 


These applications of OFDM use pilot signals and training symbols to conduct 
channel equalization. In the next chapter, we see how channel equalization can be 


conducted without the use of these methods. 
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Hil. BLIND CHANNEL ESTIMATION IN OFDM WITH ZERO 
PREFIX 


A class of blind equalization algorithms has been presented by P. P. Vidyanathan 
et al. in a number of papers [1], [5], which can easily be extended to OFDM with ZP. 
Their argument is presented in this chapter. An alternative derivation of their results is at 


the basis of further properties of this algorithm, as is presented in Chapters [V and V. 


A. BLIND EQUALIZATION USING A ZERO PREFIX 


In the method described in [1], the CP used in the standards discussed in Chapter 
II is replaced with the use of a ZP. Instead of filling the guard interval with data from the 
end of the data block or the last L samples of the data block, the data is zero-padded with 
L zeros during the guard interval. As long as the guard interval is still longer than the 


maximum time spread of the channel, there is no ISI. 


In ZP OFDM the data is transmitted as a sequence of discrete time OFDM 


symbols x[0], x{1],...,x[],..., defined as 


| x [n] 


x,[n] 


xln}=| x,[n]} €Cc™ (14) 
0 








0 
where M is the length of the FFT and IFFT, and P=M +L is the overall length of the 
OFDM symbol, including the ZP. Then each input block is in a subspace of C’™'. If 


there is no noise, the received signal 


n : 
yln]= ec” (15) 
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is also in a subspace of C’". With this in mind, it is possible to show that the channel 


impulse response c[/] can be estimated from the subspace of y[n] without knowledge of 


the input signal x[n]. 


B. CHANNEL ESTIMATION 
The general principles behind the method proposed in [1] are discussed in this 
section. Let the n-th data block be defined as 


soln] 
s[n]=| (16) 
Syl] 


The data block is then encoded with a matrix R,<¢C™”™ which in the cases of 


OFDM modulation is associated with the IFFT. 








s[n] ——>~ Ry [> 


ADD ba Ne 
ZEROS 























Figure 3. The signal s[m] is encoded with the R, matrix, then L zeros are appended. 
After [7]. 


After precoding, the transmitted blocks x[n] are defined as: 


R, 
x[n] = s[n]. (17) 
0 
Notice that x[n] can be written as 


a(n =[xp[n],x,[n],-...%y [7], 0,...0] (18) 
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with the zeros indicating the zero prefix. Although the appended zeros really represent a 


suffix, for clarity of discussion, they will continue to be referred to as the zero prefix. 


In the section that follows, it can be seen that the received signal has a particular 
structure from which the channel’s characteristics can be extracted without knowledge of 


the transmitted data. 


i Ee Received OFDM Symbol 


The output from the channel in the time domain is the convolution of the impulse 


response of the channel with x(t). Since there is no ISI due to the ZP, each block can be 


examined separately, and the received signal for the n-th block becomes 


L-1 


yin] =) clk]x,_,[7] £=0,1...,P—1 (19) 


k=0 


with c[k], k=0,...,L—1, the impulse response of the channel. Then it is easy to verify 








that 
| c[0] O = OO | 
c[1] c[0] 
meet: ge ae age 
Ae (20) 
: 0 c[L-1] -:- c[0] 
Xy— [1] 
y,aln] : 0 efi] 
| 0 0 ++ c{L-1]| 
which can be represented in matrix from as 
yln] = Ax,[n] = AR, stn] (21) 


with y[n] a Pxl vector, x,[n] , the data only version of x[n], a Mxl vector, and A a 


PxM matrix. The matrix A contains the information of the channel and it has an 


interesting property. Given any complex variable z , it can be shown that 


13 





| [0] 0 - 0 | 
c{l] c{0] 
: cll] see : (22) 
aN ne Oe ae aie. za iC) ec sateen tem Neto 
: 0 asi c{l] 
[ 0 0 ve tL -11| 





with C(z), the transfer function of the channel, defined as 


C(z) = Selle" =Z fell} . (23) 


The right hand side of (22) becomes 
CO) Ta). ame (24) 
Using the results above, we see how the transmitted data s[n] can be recovered 


from the received data yin) assuming that there is no noise. Let 
s{n]=Eyln], Eec™?, (25) 


with E,, an arbitrary M x P matrix to be determined later. Substitute (21) into (25), and 


the following equation results: 


s(n] = EAR, s{n] 


. (26) 
=[ E ]/ AJ[R,]s{n] 
Now, let E be of the form 
T 
£0 
2 ene (27) 
T 
€u-t 
where e,’ is defined as 
T 1 -l =) -(P-1) 
Co = S— ee Oe Prey Zee ,| k=0,....M-l1 

k ua: k k k ] (28) 


If a new matrix is defined as 
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C(m) Ov 0 


0 ae stg : 
| ow : (29) 
: Bg 0 
0 0 Clpy-1) 
then it can be seen that 
L. oo 
EA=A;, : 1A 
ieee 2 Fea Pate 
1 Po: parr 
= Nos Pa : : : . (30) 
I <p Cua 
= View 


Substituting the results of (30) into (26), we reach the conclusion that 
Sf] = Vy RE]. (31) 
The transmitter-receiver structure in (25)-(31) is more general than standard 


OFDM. It can be easily seen that the choice 
p.=e “ k=0,....M-1 (32) 


will yield the standard OFDM. In this case, the operation can be done using FFTs if the 


R, matrix is chosen carefully; namely, 
Revd: (33) 
2. Channel Estimation from a Block of Received Data 


If the received vectors are organized column-wise into a matrix and the number of 
received vectors is much larger than the length of the OFDM symbol such that 
Y=| yO], yf), ..., yl¥-] JP, (34) 
then it follows that Y isa PxJ matrix of the form 
Y=AR,[sO], sl], ..., s[J]]. (35) 
From (35) it can be seen that Y is not a full rank matrix since the columns of A 


in (22) span an M -dimensional subspace of C”. Therefore ,a matrix V exists such that 


15 


VY =0. (36) 
Since the matrix of transmitted data is most likely full rank, the matrix V also 


satisfies VA =O or, more explicitly, 








| [0] 0 vee 0 
c{] c[0] : 
; cf] : 
c[L-1] : te 0 
V =0. (37) 
0 c[L-1] --- — c[O] 
: 0 we c[l] 
| 0 0 ++ LI] 


The matrix V , which can be computed from the null space of Y , contains the 
information of the channel, as can be seen next. 


To extract the channel information from the V matrix, several more steps are 


required. Let G be a Px P matrix, featuring a structure similar to E in (27), 





G= : (38) 








where 


k =0,...,P—1 (39) 
P: * P; 
Although this is more general, OFDM would correspond to the choice p, = e”” for 


some w,. Now, let the G matrix be used for post-processing of the received signal. For 


the illustration of this post-processing, a new y'[n] is defined: 
16 


y[n]=GAR,s[n]}. (40) 
Notice that (38) and (40) yield a similar receiver to the one defined by (25) and 


(27), only with an extended set of frequencies at the demodulator. 


Now, substituting (38) for G into (40), we obtain 











£0 
i 
e 
y'[n]= éy, |AR sin]. (41) 
T 
eu 
T 
Ep 
If V,,.., is defined as the first M columns of G, A,, is defined in (29) and A, is 
defined 
C(ipy) O° 0 
ees 0 ie 3 og (42) 
ot “Ee, 0 
0 -- 0 C(p,,) 
then (41) can be written as 
A 0 || Vi. 
ynj=| ye R,sfn]. (43) 
_ OAL LV iw 
Since R, was chosen as V,/!,, , (43) can now be rewritten as 
[a] Sah * | sfal (44) 
nj= s[n], 
= 0 A, Br 
where 
B=V,.,R, (45) 


Recall that V is an M x P complex matrix such that the product of V and the Y 
matrix defined in (34) is equal to zero. It follows that 
V yin) =OVn. (46) 
This implies 
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ean 
V =0 (47) 
0 A,||B 


v vil |2o 48 
[V, a (48) 


or, by properly partitioning V , 


To satisfy (48), the V matrix can be chosen 
v=[V, V,]=[A,BA, -7], (49) 
where J is an LxL identity matrix. This V matrix must also be equal to zero when 


multiplied with the Y matrix, meaning that 
[ALBA;, I? =0, (50) 
Y, 


where Y, is a matrix of the first M rows of Y, and Y, is a matrix of the last L rows of 
Y . From (50) it is easy to see that 
A, BAyWY =Y,- (51) 
Using linear algebra on (51), we get 
A, BAYY, =¥,, 
A BAY Y¥o =Y%» 
A BAS AVY, HY, 4: * On 





Let 
* ok -1 
Z,=VY, (oY ) (53) 
where 
Ant Cee (54) 
Equation (52) then becomes 
Ay B= Z Ape (55) 


If the structure of (55) is examined more closely and the k -th row of both sides 


are taken, 
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Cline) Dey: 28 Dey =| aie tee al i (56) 


then \ 
CCP) | Biss tees Desay |=) ese Dads Zp C(P)s + eT LG ee a (57) 


With one last manipulation, the frequency response for each of the frequencies 


associated with the M point IFFT is 


1 b 
C(p,) a 
C( Px) Zk 
1=0,1...,.M—-1 ; (58) 
k =0,1,...,L-1 


The C(p,,,,) term is a constant associated with all frequencies of the channel 


response. Additionally, for robustness in the presence of noise, the channel response can 


be averaged over the k rows of the matrix. 


In conclusion, the channel frequency response, apart from a scaling factor, can be 


computed from a block of received OFDM symbols, as in (34), as follows: 


Step 1: Let V, 305 V, be the first Lx M and the last LxL blocks of G in 


MxM 


Equation (38). 


Step 2: Compute B = V, 


LxM 


R,. 

Step 3: Let Y, and Y, be the first M rows and last L rows of Y , respectively. 
Step 4: Compute Z, =Y,¥, ban ie 

Step 5: Call z, and b, any row k of Z, and B. 


Step 6: C(p,), C(p,), ..., C(,_,) is proportional to ratio of rows z, and b,. 


The channel frequency response was calculated using a MATLAB program that 
implemented the above method. The next chapter discusses the adaptive method used for 


real-time processing of an OFDM signal. 
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IV. RECURSIVE IMPLEMENTATION 


A. OFDM WITH CYCLIC PREFIX AND OFDM WITH ZERO PREFIX 


The blind channel equalization algorithm presented in Chapter III is not recursive, 
and a large number of received OFDM symbols are required to compute the channel 
estimate. In this chapter, a recursive method of computing the channel estimate for every 
block received using an LMS estimator is explained in detail. For the purposes of this 
portion of the thesis, a received signal similar in structure to that used by IEEE 802.1 1a is 
used. In this standard, the length of the IFFT M is 64 samples, and the length L of the 
cyclic prefix is 16, for an overall OFDM symbol length P of 80 samples [6]. 


The method of blind equalization described in Chapter III required a ZP. In this 
chapter, we will see that the CP OFDM symbol used in all standards can be reframed as a 
ZP OFDM symbol by simple preprocessing at the receiver. 


To achieve this goal, a standard CP OFDM signal, y,,[n], can be made into a ZP 
OFDM signal by subtracting a delayed version of the signal to get 

Yzpl1] = YopLn]— Yep[n—-M ]. (59) 

As shown in Figure 4, it is easy to see that the CP terms cancel since the data in 


the CP is the same as the data at the end of the OFDM symbol. To recover the old 
symbols, the new OFDM symbol and the phase shifted old symbol are simply added. 





| CP 1 | OFDM SYMBOL 1 | CP2 | OFDM SYMBOL 2 | 





CP 1 OFDM SYMBOL 1 CP 2 


OFDM SYMBOL 1 WITH 
PHASE SHIFT zp NEW OFDM SYMBOL 


Figure 4. | Adelay of M points results in a phase shift of the first OFDM symbol with 
a ZP and new symbols subsequent to that. 
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Now the channel can be estimated using the algorithm discussed in the previous 
chapter and applied to the received OFDM symbol. The rest of this section describes the 
LMS filter method used as if a zero prefix OFDM symbol had been transmitted vice the 
CP OFDM symbol that was actually simulated. 


Initially, an LMS estimator was derived directly from the off-line method 
described in Chapter III. This method is not robust, and the theory behind it is not 
presented in this thesis; however, some of the simulation results have been included for 
comparison purposes. In what follows we present a different approach to blind 
equalization for ZP OFDM, discussed in [9], which is more suitable to recursive 


implementation. 


B. PROPOSED APPROACH 


In this section, a different approach to blind equalization of ZP OFDM is 
investigated based on standard properties of the DFT. The results of this section are 


based on [9]. 


The purpose of this approach is an improved estimation method that can be 
computed recursively with an LMS. The resulting method is shown to be effective in 
cases where the channel is initially poorly estimated or slowly changing with time. It 


also had an added benefit that is discussed in the next section of this chapter. 


In order to introduce this approach, define a vector x of length 2M such that 
x[n]=0, forn=M,...,2M —-1 (60) 


and a second vector w of the same length 2M of the form 


0 if0<n<M-1 
wen] = ; : (61) 
1 iffM<n<2M-1 
Clearly, no matter what x[n], n =0,...,M@ —1 is, we have 
x[n]w[n]=0 for alln=0,...,.2M -1. (62) 


Since multiplication of two finite length sequences corresponds to the circular 
convolution of the respective DFTs, we can write (62) in terms of W[k] = DF T {w[n]} 
and X[k]= DFT {x[n]}. 
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The DFT of the w vector is defined as 
Wikl= Sie 2m", (63) 


which is a geometric series. If /=n-—M >n=I1+M,, then (63) can be rewritten as 


N-l _ 22,4, —(—1)* 
Se toes M) Fer ea bs GU, for k =1,...,2M —1 

WIk] =3 t= fae ' es 
M, fork =0 


It is easy to see that W[k]=0 when k is even, and in particular 


if k odd 
-j—k 
e N 


W[k]= if k #0, even. (65) 


(< 
0 
M ifk =0 


Now, circularly convolve the DFTs of the two vectors x and w to get 
M-| 
> X(Wik-1]=0, (66) 
1=0 
where the shifted represented by k—/ is a circular shift. If k is chosen to be any odd 


frequency, it can be shown that W[k —/] is not equal to zero only whenk =/ or / is even. 


Using this conclusion, we see that when k is odd (66) becomes 


3 X(IJW[k —1]+ X[k]W[0]=0. (67) 


1=0 
leven 


If the value for W[0] is substituted into (67) and the equation is rearranged, then 


MX{k] ay, Wik -J]X[1]. (68) 


1=0 
leven 


The relationship between the DFT of the transmitted signal and the DFT of the 
received signal is just the inverse of the channel frequency response; that is, 


Y[k] 
X[k]=——. 69 
[k] Clk] (69) 


For any odd value of k , (69) can then be substituted into (68) to get 
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1 = 1 
——=-| ) Wik -“Y— |. 
CIk] » oe aro 


leven 


MY[k] (70) 
If (70) is further deconstructed, with the / values associated with the data on the 
left-hand side, the / values associated with the pilots on the right-hand side, and an 


arbitrary odd value of k, then 


M-1 1 1 reat 
W[k -1]Y(]——+ MY[k] ——=-| Y WIk-11X,,,,.{/] |. 71 
py ie ont wars 2 [K —FVX ys LE] (71) 


leven leven 


With a few terms further defined, (71) lends itself to an LMS-type 


implementation; that is, 





@.,=[WIKIVIO], WI(k-2)1¥I2], ..., MYLK]] 
0 aay ar ae a ] (72) 
Clo] C12] CIM -2]' CIk] 


M-1 
Vex SS » W[k —1X sioislt] 


1=0 


Now, from [9], the LMS estimator can be written as 
ae = 8 ua 25 up. (a ., = vz) : (73) 


Cc: SIMULATION RESULTS 


Each of the methods of implementation for the blind estimation of the channel 
were run using 100 trials at each data point, regardless of the type of test being run. 
Additionally, the number of OFDM symbols transmitted was set at 10,000. Using IEEE 
Standard 802.11a format, we have 52 BPSK or QPSK symbols per OFDM symbol, for a 


total of 52 million symbols per data point. 


The number of pilot signals was reduced from four pilot symbols per OFDM 
symbol to one pilot symbol per five OFDM symbols. The purpose of transmitting the 
pilot symbol was to determine the channel constant for all three methods of 
implementation. However, the number of pilot signals sent could be further reduced 


depending on channel conditions. 
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The values of wz, the LMS algorithm step size, were determined experimentally to 
ensure convergence of the adaptive algorithm. For the initial LMS implementation the 
step size 42 was set to a small value (107’), and for the second DFT-based LMS 
implementation, jz was set to a larger value (3x10~). These values can be adjusted to 


ensure the algorithm converges or to adjust the speed at which the algorithm converges. 


1. Performance in the Presence of Multipath and AWGN 


To test the performance of the proposed DFT-based adaptive channel estimation 
method and the adaptive method based on the theory discussed in Chapter III, the 
respective performances were compared using the initial off-line, or block, blind 
estimation method. All three methods were tested with AWGN. The tests were run with 
SNRs from zero to 50 dB. Algorithm performance was determined by calculating the 
symbol error rate (SER), in errors per symbol. The results for the DFT based adaptive 
algorithm with an LMS step size of 4 =107 and the block method is shown in Figure 5. 


The channels are modeled as random complex channels with an impulse response of the 


same length as the guard interval ( ZL =16 in this case). 
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Figure 5. | QPSK SER versus SNR with a step size of w=10"'. 


For this thesis the SNR is defined as 


_ Es*()] 
E[n*(f)] 





this factor is 0.692 [6]. 


(74) 


where s(t) is the signal and n(t) it the noise. It is related to E,/N, by a factor of B/R, 


where B is the bandwidth of the channel, and R, is the data rate. For 802.11a QPSK, 


It can be seen that, in this example, the proposed LMS method has better 
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performance than the block processing approach presented in Chapter HI for SNR less 
than 18 dB. At 5 dB, the proposed method out performs the block method by 
2.423x10" errors/symbol, the maximum delta between the two methods. This is likely 


due to the fact that the block processing approach estimate is averaged over the k rows 


of the matrix, as discussed in the previous chapter. Specifically, there are k estimates of 
the channel, one derived from every row of the matrix, and these estimates are averaged 


together before being applied to the received signal. 


The results for the adaptive method discussed in Chapter III are not presented 
since the algorithm never fully converged to the correct channel. This is due to the fact 
that the matrices that the LMS method is based on are ill-conditioned. The block method 


uses singular value decomposition to overcome this limitation. 
The next simulation was done under the same conditions except the value of the 
step size was changed to 44=3x10° to determine the effect on the proposed algorithm. 


The results of this change can be seen in Figure 6. The proposed adaptive method works 


poorly compared to the block estimation method. 
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Figure 6. QPSK SER versus SNR with a step size of uw =3x10°. 
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25 Non-Stationary Channel Performance 
The next test of the proposed implementation was to examine the effect that a 
non-stationary channel has on its performance. The non-stationary channel was 
simulated using a basic random walk. The channel was initially a set of four complex 
impulses distributed within the length of the guard interval. Every OFDM symbol, this 
channel was updated to a new value by a random walk. The random walk step size has 
been calculated as a percent of the channel power. The channel coefficients are obtained 
iteratively from 
C.,, =C, + BELC, JAC (75) 
where AC is a random vector consisting of [1,0,-1] , # is ascaling factor, and E[C?] is 


the channel power. 


The block method of channel estimation was expected to perform poorly, based 
on the fact that it assumes a stationary channel over the entire data set of 10000 received 
blocks. With the small LMS step size =10', the adaptive method performed 
comparable to the block method. The small step size affected its ability to converge on 
the actual channel estimate as each OFDM symbol was received. 

To ensure that the poor performance was due to the LMS step size, the simulation 
was also run using “4=3x10~. The results of this simulation are shown in Figure 7, 


where SNR is defined in Section 1 of this chapter. 
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Figure 7. |. QPSK SER versus random walk step size with an LMS step size of 
u=3x10>. 


As illustrated in Figure 7, the LMS method based on the blind equalization 
algorithm in Chapter III is very sensitive to the LMS step size. The proposed DFT based 
adaptive method performs better than the block method and the Chapter III LMS method 
fails to converge to a solution. The step size for the proposed DFT based adaptive 
method can be adjusted to give faster convergence on a non-stationary channel; however, 
the gains are marginal, since the block method reaches 10% SER for a random walk step 
size of 0.21 percent of channel power, while the alternative LMS method reaches 10% 
SER at 0.23 percent. This can also be seen in Figure 8, where the alternative LMS 
method and block method have been plotted using a logarithmic scale. The error rate for 
the block method is better than that of the alternative LMS method until a random walk 


step size of 0.08 percent of the total channel power is reached. Although there are several 
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data points at which the block method outperforms the alternative LMS method above 


0.08 percent, in general, the alternative LMS performs better. 
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Figure 8. | A comparison of the block and proposed adaptive method of calculating the 
channel using a step size of 44 =3x10~. These results were generated 
using 500 trials at each data point. 


These results also fail to take into account the fact that the larger step size affects 


the performance in AWGN, as seen in Section 1. 


3. Effects of Errors in the Initial Channel Estimate 


The simulations performed in the preceding sections assumed that there was a 
perfect initial estimate of the channel conditions based on the reception of a preamble. 
The performances of the algorithms based on a poor initial channel estimate are 
investigated in this section. Only the two LMS-based methods are tested in this section 


since the block method does not rely on an initial channel estimate. The error in the 
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initial estimate was simulated by taking the original channel estimate and multiplying 
each frequency of the channel by a different random real number. The initial estimate is 
given by: 

c 


init, j 


= Cre iV (76) 


act,j j 


where C. 


init 


is the channel estimate given to the simulation, C,,, is the actual channel 
impulse response, and v, ~ N(1,r) with r defined as the seed channel variance. The 


performance metric for this algorithm was based on the number of OFDM symbols it 
took for the algorithm to converge to less than a 10% SER. This was calculated by 
running the algorithm until an OFDM symbol had five or less errors in it. The number of 
symbols was still set to 10,000, and if the algorithm failed to converge by then, it was 


assumed that it would not converge. The results of this simulation are seen in Figure 9. 
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Figure 9. | Number of OFDM symbols required to converge versus seed channel 
variance for a step size of 4 =10°'. 
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Both algorithms perform poorly at the lower step size of 4=10', taking over 
8,000 blocks to converge at a seed channel variance of 1.0 and failing to converge 


beyond a seed channel variance of 1.4. Assuming a OFDM symbol interval of 4.0 ws, 
based on IEEE Standard 802.11a, we find that this translates to a convergence time of 32 


ms. 


The step size was then increased to “w= 3x10 to see what effect, if any, it had 


on the convergence of the two algorithms. It can be seen in Figure 9 that the alternative 


LMS method converges much faster than the original LMS method. 
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Figure 10. Number of OFDM symbols required to converge versus seed channel 
variance for a step size of 4=3x10°. 
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With the larger step size, the alternative LMS method converges between 20 and 
24 ms for seed channel variances of 1.1 to 2.0. The original LMS method shows slightly 
degraded performance at the larger step size. Due to the fact that the algorithms do not 


always converge, the plot is not smooth through seed channel variances of 1.3 to 2.0. 


This chapter investigated a DFT-based adaptive method for channel estimation. 
The next chapter discusses an algorithm for the estimation of data in null subchannels 


that is also based on the DFT. 
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V. DATA RECOVERY FROM A FADED SUBCARRIER 


A. PROPOSED NULL ESTIMATION METHOD 


One of the issues in OFDM is the loss of data in subcarriers corresponding to 
nulls in the frequency response of the channel. This is particularly significant in 
environments with multipath, which adds destructively at certain frequencies. The 


method presented in this chapter is based on [10]. 


The argument of this section is that, using OFDM with ZP and an FFT size twice 
the OFDM symbol length, the even and odd subcarriers of the OFDM symbol offer 
sufficient redundancy so that the lost data can be recovered. 

In order to introduce this concept, let the transmitted OFDM symbol be 

x =[x{0], x{1],..., x1 -1],0,...,0] (77) 
where x is a vector of length P, which includes the ZP of L zeros at the end. Now 
define the M point DFT and the 2M point (zero padded) DFT as 

X[k] = DFT {[ x{0],....x1M —I]]}, k =0,....M—-1 


, (78) 
X yl] = DFT {[x{0],...,x[M —1]],0,0,...,0}, m=0,...,2M -1 
where M zeros are added to x before the 2M point DFT is taken. 
Then, as discussed in (66), for any odd k we can write 
2M -1 
> X,yJWIk —1]+ X,,,[k]W[0] =0. (79) 


1=0 
leven 


Since ¢ is even and k is odd we can let k =2m-+1 and rewrite this equation as 


Y2m+l] z 
M EST 2, W[2(m—-1) +1]X,, [1]. (80) 


Now, set an arbitrary threshold and use it to partition the transmitted data into two 
parts such that: 
X[ =X, [1+ X_[I] (81) 


with 
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X(]_ if |Cl]|> Thresh 
X [I= 
0 if |C[l]|<Thresh 


Xl] if |C[l]|< Thresh 
xX _[l]= 
0 if |Cll]|>Thresh 


(82) 


where Thresh is an arbitrary threshold. 


In other words, X_[/] indicates the subcarriers received at low power (at or close 


to a null in the channel) and X [4] all others. 


Now, if the partitioned data is inserted into (80), it becomes 


Y2m+l]_ 1% ; ie ; 
Cons le ee * AB 
Define 


j 


w)ln] = IDFT {W [2k +1] =-e 7H" 


(84) 
x,[n] = IDFT {X,, [1]} 
Using the above definitions in (83), we get 
Y[2k +1] _ 
Cl2k+ i} DFT {wy[n]x,[n]} + DFT {wo[n]x_[n]} . (85) 


If the term containing x_[n] on the right hand side of the equation is isolated, then 


Y[2k +1] 
DFT {wo[n]x_[n]} = apr eel ie DFT {woln]x,[n]} . (86) 
If the inverse discrete Fourier transform (IDFT) of both sides of (86) is taken, 
then 
Y[2k +1] 
= [DFT , ———— > - : 87 
woln]x_[n] Cl2k + | woln|x,[n] (87) 


Now, dividing both sides by the w,[n] term to solve for the terms of the 


transmitted signal below the threshold, we get 


I Y[2k +1] 
S558 gp) er. 
x[7] win] ae x,[7] 





and, finally, 
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X_[¢|=DFT {x _[n]}, ¢=0,...M—-1 (89) 


is the data recovered at, or close to, the nulls of the channel. 


B. SIMULATION RESULTS 


To test the null estimation method described in Section A of this chapter, a test 
similar to those in Chapter IV were run, with several significant differences. Instead of 
using the OFDM symbol format described in IEEE Standard 802.11a, no nulls or pilots 
were included in the OFDM symbol. This was done in order to test the proposed concept 
of null estimation without having to worry about the effects of the nulls in the OFDM 


symbol. Additionally, knowledge of the channel was assumed. 


The null estimation algorithm continued to use the delay of one block to achieve 
the effects of a ZP while transmitting the symbol with a CP. To test the algorithm, 


transmission was simulated over a noisy channel, where the noise was AWGN. 


For testing, we used a family of random channels all with zero frequency response 


at a specific frequency @,. An example channel response is shown in Figure 11. 
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(a) 


w (Radians) 
One of a family of channel frequency responses, all having nulls at 


Figure 11. 
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Figure 12. A comparison of the QPSK symbol error rate between a normal OFDM 
receiver algorithm and a null estimating OFDM receiver algorithm. 


The performance of the standard algorithm (solid line) is compared with that of 
the proposed algorithm (dashed line) in terms of symbol error versus SNR in Figure 12. 
The two algorithms have similar performance up to 25 dB SNR, although the null 
estimating OFDM receiver algorithm performs slightly better overall. Above 25 dB 
SNR, there is a sharp drop in the error rate from the null estimation algorithm, while the 
normal OFDM receiver algorithm remains constant at 7.8x10~ errors per symbol due to 


the information lost in the deep fading sub-carriers. 


This chapter presented an algorithm for the estimation of the data lost in null 
subchannels and simulation results from the implementation of this algorithm. Chapter 
VI provides a summary of the work presented in this thesis, significant results, and 


recommendation for future work. 
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VI. CONCLUSION 


The goal of this thesis was to investigate an adaptive implementation of the blind 
channel equalization method proposed in [1]. This adaptive method was implemented 
using the current OFDM symbol architecture from IEEE standard 802.1la. Once the 
adaptive method based on the work in [1] was simulated, an alternative method was 
developed with more robust performance characteristics. The algorithm for the second 
method led to the development of a proposed method of null channel estimation, which 


was also simulated in this thesis. 


A. SUMMARY 


This thesis presented a subspace-based method of blind channel estimation of 
OFDM modulated signals that used the zero prefix as the guard band interval. Based on 
this, a method of adaptive blind channel estimation has been proposed. Furthermore, it is 
shown that the redundancy in ZP OFDM can be exploited to recover data from deeply 
faded subcarriers. The blind channel estimation methods were implemented in MATLAB 
and then compared for performance for several different channel conditions. Both a 
noisy channel with AWGN and a non-stationary channel were used to compare the 
results of the three methods of channel estimation. Additionally, two adaptive methods 
were tested in the presence of errors in the initial channel estimate to see if they could 


recover the original channel estimate. 


Finally, the algorithm for data recovery from faded subcarriers was coded in 
MATLAB and simulated. Since in this case we assume full knowledge of the channel, it 
was only tested in a channel with AWGN. The results of this simulation were compared 


to those of an OFDM system that did not perform null estimation. 


B. SIGNIFICANT RESULTS 


The method of blind channel estimation proposed in [1] can be implemented in an 
adaptive form. However, this method is very sensitive to the step size used in the LMS 


algorithm. The algorithm did not converge in a channel with AWGN with a relatively 


4] 


large step size (44=3x10~ in the example). Similarly, the algorithm performed poorly 
with this step size for both the non-stationary channel simulation and the incorrect initial 


channel estimate simulation. 
The proposed adaptive method is more robust than the original. For the AWGN 
noisy channel with an LMS step size of 44=10~’, the system converged to a lower error 


rate in simulated low SNR environments than the off-line, block method. This can be 


seen in Table 1. 


Table 1. Comparison of calculated SER at various SNR for the blind estimation 
methods using a step size of w=10". 























Proposed DFT 
Block Method Based Adaptive 
enh) (Errors/Symbol) Meied 
(Errors/Symbol) 
0 5.799 x107! 4.543107! 
A 4.513x10" 2.09010" 
10 1.741x107! 6.203 x 10 
IS 3.895 x 107 2.155x107 
20 2.640107 5.982 x10° 
25 1.390x10° 4.837x10~ 

















With the larger step size, the block method outperforms the adaptive LMS 
algorithm. However, there are several advantages to the larger LMS step size, which can 
be seen in the results from the other two simulations. The alternative adaptive method 
with the larger step size outperforms the block algorithm during the non-stationary 
channel simulation above a random walk step size of 0.08 percent of channel power by an 
average of 0.0175 symbol errors/symbol. For the bad initial channel estimation 
simulation, the larger LMS step size outperforms the smaller step size significantly, 


which can be seen in Table 2. 
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Table 2. Comparison of convergence speed for various seed channel variances for 
the original and alternative adaptive methods with an LMS step size of =3x10°. 

















Sat Channel Adaptive Method Alternative Adaptive Method Time 
Wanance (Number of OFDM (Number of OFDM Difference 

Symbols) Symbols) (ms) 

0 1 1 0 

0.5 1001 4.25 4 
1.0 9800 4009 ike 

3 Fails to Converge 5049 00 

2.0 Fails to Converge 6031 00 














Perhaps the most significant result is the fact that the proposed adaptive method 


can be tuned for performance based on the channel conditions that may exist. 


The recovery of data transmitted over deep fading channels was shown to be 
possible using the algorithm described in Chapter V. For relatively high SNR (25 dB in 
this example), the proposed algorithm clearly outperformed the standard receiver which 


cannot recover the faded data. 


C. RECOMMENDATIONS FOR FUTURE WORK 


There are numerous areas available for future work with this thesis. First, the 
proposed algorithm should be evaluated using actual channel characteristics to determine 


if it is viable for use in an actual communication method. 


Secondly, the adaptive method should be tuned to perform based on channel 
characteristics. In a normal OFDM modulation scheme, such as IEEE Standard 802.1 1a, 


a known preamble is sent before data transmission begins. The channel characteristics 
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determined from the preamble should then be used to determine the desired LMS step 
size. By modifying the step size, the performance of the algorithm can be enhanced for 


noisy or non-stationary channels. 


Third, the null estimation algorithm should be tested using an OFDM symbol 
from a current standard and integrated with channel estimation method. If knowledge of 
the channel were known before transmission, no data would be sent over channels with 
deep fading. Once a channel has been determined to have sub-carriers with deep fading 
characteristics, null estimation should be performed to determine if the algorithm will 


perform as well with only an estimation of the channel characteristics. 
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